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I. INTRODUCTION 



There appear to exist some misunderstandings concerning the determination by our 
group of the pion-nucleon coupling constants in multienergy partial- wave analyses of pp [|TJ , 
combined pp and np 0, and pp [g] scattering data. We have summarized the recent deter- 
minations of the pion-nucleon coupling constants in Table |[ 

In their latest analysis of 7r ± p scattering data the VPI&SU group has obtained a value 
for the pion-nucleon coupling constant consistent with our values M . Since the values f^ Nn 



0.075 found by us (see also Refs. P,|],|7|,|8||l ) are consistently lower than the value f% N% 



0.079 found in the Karlsruhe-Helsinki analyses of u p scattering data [10.11,12], a number 



of papers fl^JT^JT^Jl^Jl7|Ji^JT9|] have been published commenting on our results. In several 
of these papers one attempts to find an explanation for the difference by trying to point out 
alleged shortcomings and so-called systematic errors pertaining to our method of analysis. 
In view of the existing confusion, we have looked again into these and other matters and 
carefully reexamined possible sources of systematic errors. We felt it would be useful to 
supply some additional information to further clarify these issues and to address a few 
points where some of these papers criticizing us are in error. 

A very important point that has not been recognized and appreciated enough is that 
the basis of our accurate determination of the pion-nucleon coupling constants and of the 
nucleon-nucleon phase shifts is our multienergy partial-wave analysis of all nucleon-nucleon 
scattering data below Ti a b = 350 MeV. This multienergy partial-wave analysis is much 
more sophisticated and computer intensive than other similar analyses of nucleon-nucleon 
scattering data. The energy dependence of the phase shifts in our analysis is described in 
a much better way. As a consequence, we end up with a multienergy partial-wave solution 
which gives an excellent fit to all nucleon-nucleon scattering data below Ti a b = 350 MeV. The 
phase shifts and pion-nucleon coupling constants as determined by this multienergy partial- 
wave analysis are the best possible values for these quantities that can be obtained from 
these data. We strongly feel that statements about the pion-nucleon coupling constant made 
with the help of single-energy analyses or, even worse, analyses of individual experiments 
with the help of potential models cannot, in any way, be compared in quality with the 
results of these multienergy analyses. The reason is that in such analyses at one particular 
energy the information about the energy dependence of the phase shifts due to one-pion 
exchange cannot be incorporated. In a multienergy analysis the complete database is used 
and energy-dependent constraints are properly accounted for. This results in a much better 
and more precise value for the coupling constant, with a minimal model dependence. For 
a proper determination the whole database is needed and arguments based on studies of 
specific experiments are not so reliable. 

In our discussion, we take mostly the Nijmegen multienergy analysis of all pp scatter- 
ing data below Ti a b = 350 MeV [[TJ as specific example. The value for the ppn° coupling 
constant obtained in this analysis is undoubtedly the most compelling evidence for a low 
pion-nucleon coupling constant, since in this case the required theoretical input is rather 
small and practically model independent. Therefore, the possible systematic error in this 
case is very small. Actually, the statistical error on the npir^ 1 coupling constant obtained in 
the np analysis is smaller than the statistical error on the ppn coupling constant, although 
the np data are less accurate and less varied than the pp data. The reason, we think, is 
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that in np scattering pion exchange can be probed more easily, since in pp scattering one 
always has to deal with the infinite-range repulsive Coulomb potential and one first has to 
strip the pp data from this long-range Coulomb contribution in order to get a handle on 
pion exchange. Nevertheless, in the pp analysis we feel almost certain that the systematic 
error on the ppn coupling constant is small, since in this case we explicitly investigated 
all thinkable sources of systematic errors and found no significant effects. We stress that 
our emphasis on the ppir coupling constant does not mean that we have doubts about the 
correctness of the determination of the pion-nucleon coupling constants in the np and pp 
partial-wave analyses. It simply is more difficult to do a similar thorough study for the np 
and pp cases, where the amount of theoretical input into the analyses is larger. For instance, 
in these cases it is necessary to make from the start some theoretical assumptions about the 
validity of charge independence. 

In the following sections of this paper we will report on our search for systematic errors 
in the value for the ppir coupling constant in the pp analysis. Part of the reason we have 
been able to do this comprehensive investigation is that recently a lot of computing power 
has become available to us. We begin the next section by discussing some statistics relevant 
to the analyses. Then we will subsequently discuss the influence of form-factor effects, the 
sensitivity of the different types of observables to the pion-nucleon coupling constant, and 
its determination from individual partial waves and from data in different energy ranges. 
Next, we will investigate more closely some particular pp and np scattering experiments that 
are brought up in connection with the determination of the coupling constant. Finally, we 
spend a few words on the Nijmegen analysis of pp scattering data. 



II. STATISTICAL CONSIDERATIONS 

The Nijmegen database of pp scattering data below T lab = 350 MeV contains at present 
Nobs — 1656 scattering observables. Since there are 119 experimental groups with a finite 
overall normalization error as well as 12 additional normalization parameters (from angle- 
dependent normalizations) we have a total of iV^at = 1787 pp data. In the Nijmegen partial- 
wave analysis these data are fitted with iVp ar = 22 model parameters, which includes the 
ppir coupling constant. Because there are also 22 groups with a floated normalization 
the number of degrees of freedom is N df = 1612. If the database is a correct statistical 
ensemble and if the theoretical model is correct, one expects (x m in) = -^df =t v^Vdf = 
1612 ± 57. In our latest analysis we reach x m m = 1786, which is only 174, or 3 standard 
deviations, higher than the expectation value. This difference is at least partially due to 
small theoretical shortcomings in our model. This implies, therefore, that there is still some 
room for theoretical improvements in the pp partial-wave analysis. We have investigated the 
statistical quality of the final pp data set by calculating the momenta of the theoretically 
expected ^-distribution and comparing them to the ones actually found in the analysis. The 
results are presented in Table JJ. For details about the statistical tools used in the Nijmegen 



partial- wave analyses we refer to Ref. ||20|| . The second and higher central moments found 
in the analysis are in excellent agreement with their expectation values. This shows that 
the statistical quality of the pp data set used is very good. It shows that our analysis with 
its Xmin/-^df = 1.108 is already quite good. Nevertheless, we hope that a significant drop in 
Xmin can still be obtained giving better agreement between and xLin- m this way our 
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partial-wave analysis becomes a tool, because it can decide on the quality of the proposed 
improvements in theory. 

We have demonstrated that our multienergy partial-wave solution is essentially correct 
statistically This has rather strong consequences: it means that the values for the phase 
shifts and coupling constants as well as the statistical errors on these quantities as determined 
in the multienergy analyses are essentially correct. Here the statistical error on a particular 
quantity (e.g., a phase shift) is the error as obtained in the standard way via the x 2_r i se - 
by-one rule, as discussed for example in Sec. V A 2 of Ref . p0| . Including new experiments 
in the multienergy analysis will change the phase shifts and coupling constants not more 
than 1 or 2 multienergy standard deviations. The same is not necessarily true for quantities 
determined in a single-energy analysis. For example, in a single-energy np analysis around 
100 MeV there are no np spin-correlation data to pin down the e\ mixing parameter, whereas 
in the multienergy analysis the presence of spin-correlation data at the adjoining energies 
near 50 and 150 MeV also allows for an accurate determination of the E\ mixing parameter 
at 100 MeV. Therefore, the multienergy values and errors for the phase shifts and coupling 
constants are much more realistic than the single-energy determinations. 

In our latest pp analysis |J we used a parametrization for the 1 Sq partial wave different 
from the one used in previous analyses . This parametrization has less parameters but 
gives a somewhat higher Xmm- We feel, however, that it is a better parametrization, since 
the 1 S'o wave was probably over-parametrized in the older analyses. This means that results 
presented for this analysis are not necessarily also true for the older analyses, especially with 
respect to statements about the 1 So wave. In spite of a large amount of effort on our side, 
the description of this partial wave is still not as good as we would like. 

In our pp analysis we determine the ppn coupling constant at the pion pole. We find 
now /p p7r o = 0.0750(5), where the error is purely statistical. If we fix this coupling at the 
old value 0.079 the result is xLin = 1842, which is Ax^ in = 56, or 7.5 standard deviations, 
higher than the minimum. We stress again that no particular data set is responsible for 
the specific value of the coupling constant, but that all data when fitted in a multienergy 
partial-wave analysis contribute to this value. From these numbers one can get a feeling 
about how unreliable it must be to make statements on the pion-nucleon coupling constant 
with the help of potential models which fit the data with Xmin/^dat — 2, instead of drawing 
conclusions based on multienergy analyses with Xmin/^dat — 1- Nevertheless, in some recent 
papers fli7l , |l~9"| conclusions about the pion-nucleon coupling constant are drawn with the 
help of potential models that have Xmm/-^dat — 2 or Xmm — 3600 which is about 2000, 
or 35 standard deviations, and not just 174, or 3 standard deviations, higher than the 
expectation value (Xmm) = 1612. The point we want to make here is that these potential 
models can still be improved in so many different places in so many different ways, that it 
is very presumptuous to try to make any conclusions about a Ax^ in = 56 effect. All such 
conclusions are inevitably very model dependent, which seems to be generally overlooked by 
the authors criticizing us. 



III. FORM FACTORS 

Especially persistent is the suggestion of systematic errors due to form-factor ef- 
fects |fr3] , ^4] , P^| , P^| , P^| , p!U| , although we have repeatedly stressed that we determine 



4 



the pion-nucleon coupling constant at the pion pole, and that therefore form factors are 
irrelevant to the issue. We did explicitly check this by adding an exponential form factor to 
the pion-exchange potential in our analysis. We used 

F(k 2 )=exp[-(k 2 + m 2 )/A 2 ] , (1) 

normalized such that at the pion pole F(— m 2 ) = 1. For values of the cutoff mass A as low 
as 500 MeV we have found no significant changes. This can be seen in Table |TJ, where the 
results for the pp7i° coupling constant and the corresponding values for xLm are presented 
as a function of the cutoff mass Ao in MeV. Evidently there are no significant changes for 
realistic values of the cutoff mass. In spite of our statements, however, it seems to have been 



suggested recently in the panel discussion at the Adelaide conference [0 that the value of 
the coupling constant may depend critically on the shape of the form factor and that we, 
like other groups, should use a form factor of the Feynman type 



F(k 2 



A^-m 2 
A 2 + k 2 



(2) 



again normalized such that at the pion pole F(— m 2 ) = 1. The square appears because 
one takes a monopole form factor at each vertex. Let us start by saying that we strongly 
feel that an exponential form factor is more physical than a Feynman-type form factor. The 
exponential form factor follows quite naturally from Regge-pole theory and from constituent- 
quark models with harmonic-oscillator wave functions. A Feynman-type form factor, on the 
other hand, which is essentially a phenomenological regulator, has an annoying singularity 
at k 2 = —A 2 in the unphysical region. This Feynman-type form factor is chosen mainly for 
reasons of convenience. 

To our mind, it is hard to understand how the type of form factor can matter once it 
has been demonstrated, using one particular type, that the value of the coupling constant 
is determined at the pole. When comparing form factors, it should always be kept in mind 
that the form factor is related to the size of the nucleon. An exponential form factor and 
a Feynman-type form factor give approximately the same nucleon size when Ao = A2/V2. 
Reliable values of the cutoff mass are hard to determine in AW scattering since these depend, 
for instance, on which heavy mesons are included in the model [|19|| . But, in general, if a 
cutoff mass is used which affects the pion-exchange potential drastically at distances where 
the nucleons do not overlap anymore, the obvious conclusion should be that this is an 



unrealistic low value for the cutoff mass. As an example, look at Figure 4 in Ref . []rj| , where 
one is willing to accept a 20% reduction of the tensor force due to pion exchange at a distance 
of 2.5 fm. In the Nijmegen soft-core potential |H[ the exponential cutoff mass is 965 MeV 



and in the Bonn potential p2| the Feynman-type cutoff mass for the pion-exchange potential 



is 1300 MeV. In fact, it has, until very recently ||19|| , always been claimed by the Bonn group 



that a satisfactory fit to the AW scattering data is impossible with a lower cutoff mass. If 
one looks at potential models, a value of 500 MeV for a cutoff mass in a form factor is quite 
low. 

To meet all criticism, however, and to avoid new misconceptions, we have again explicitly 
checked our conjectures by repeating the pp analysis using a pion-exchange potential with 
Feynman-type form factor. As we expected, our findings were entirely similar to those 
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with an exponential form factor. We conclude therefore that neither the shape of the form 
factor nor the value of its cutoff mass (as long as it is not unreasonably low) has a significant 
influence on our determination of the pion-nucleon coupling constant. The obvious reason for 
this nice feature is to be found in the specific method of analysis which allows the extraction 
of the coupling constant from the asymptotic behavior of the one-pion-exchange potential 
in configuration space and its determination is not sensitive to short-range modifications. 



IV. DETERMINATION FROM 
DIFFERENT OBSERVABLES, 
PARTIAL WAVES, AND ENERGY RANGES 



It is an interesting exercise to investigate which particular types of observables are the 
most sensitive to variations in the coupling constant. We have repeated the pp analysis for 4 
different values of the ppir coupling constant and the np analysis, with the NNtt coupling 
constants fixed at 0.075, for 4 values of the npir^ coupling constant. The resulting values 
of 

Xmin are tabulated in Table [TV] for the different types of pp scattering observables and 
in Table [V] for the np observables. The value found for the charged-pion coupling constant 
in the np analysis is / c 2 = 0.0748(3). It can be seen that in the np analysis the difference 
between f 2 c = 0.075 and 0.079 is A^ n = 255, so this is a difference of 16 standard deviations! 
We did already mention above that the corresponding difference in the pp analysis, for the 
ppn° coupling constant, amounts to Ax^ in = 56, or 7.5 standard deviations. Furthermore, 
it can be seen from these Tables that no particular type of observable is solely responsible 
for the low value of the coupling constant, but that essentially all types oi pp as well as np 
scattering data favor a low pion-nucleon coupling constant. 

In order to investigate in what way the different partial waves are sensitive to the ppn 
coupling constant, we introduced first of all two different couplings, one for the singlet waves 
and one for the triplet waves. We then find Xmin = 1786 for 1611 degrees of freedom. The 
resulting coupling constants are 0.0753(7) for the singlet waves and 0.0750(6) for the triplet 
waves. This shows that both singlet and triplet waves favor a low coupling constant. 

Next, we introduced in turn a different ppn coupling constant for a specific partial wave 
and one for all other partial waves. This exercise was done for all parametrized waves in the 
analysis: ^o, l D 2 , 3 P , 3 Pi, 3 Pt- 3 F 2 , 3 F 3 , and 3 F 4 - 3 F 4 . The results for the coupling 
constants are shown in Table IVU For all cases excepting the 1 S'o wave the results favor 



a low coupling constant. In the case that a separate ppir coupling was introduced in the 
1 S'o channel an improvement of Ax^ in = 6.6 was found for a high coupling constant. This 
is presumably a consequence of our new parametrization of the 1 Sq partial wave which is 
evidently not perfect. When the coupling constant in the 1 So wave is fixed at 0.075 the 
results for all remaining partial waves are nicely consistent with a ppir coupling constant 
somewhat lower than 0.075. Apparently the 1 S wave enhances this to fp pn0 = 0.0750(5). 
This is probably also the reason that of the different types of observables only the differential 
cross sections favor a coupling somewhat larger than 0.075 (see Table [IV]), since differential 
cross sections are more sensitive to this wave than spin-dependent observables. There is 
thus some indication for a small systematic error pertaining to the 1 S'o partial wave in the 
pp analysis. Further study is required to find out to what extent this is to be attributed to 
friction between the data or to a flaw in our theoretical treatment of the 1 5'o channel. In view 
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of the above findings, it may be better to fix the coupling constant in the 1 So channel at a 
value of 0.075 and determine f ppn o from the remaining partial waves. We then find for 1612 
degrees of freedom xLin = 1787 and the coupling constant becomes f ppn0 = 0.0746(6). This 
is within one standard deviation from the value quoted above f pp7r o = 0.0750(5) determined 
from all partial waves including the 1 Sq. 

We can also bypass the apparent friction in the 1 Sq partial wave by removing from 
the database the data taken at very low energies below 1 MeV, namely the Los Alamos 
cross sections around the interference minimum measured by Brolley et al. |23[ and the 
Zurich cross sections from Thomann et al. The reason, of course, is that at these low 
energies the 1 5'o phase shift is very accurately known and gives a very strong constraint on 
the parametrization of the l So phase shift. The results from this 3-350 MeV partial- wave 
analysis are presented in Table |V1J| . As expected, we find a somewhat smaller coupling 
constant f p = 0.0743(6) instead of 0.0750(5). If we fix the coupling constant in the 1 S 
channel at 0.075, we find in the 3-350 MeV analysis f pp7T o = 0.0744(6). So there are strong 
indications that the value for the neutral-pion coupling constant as determined in the pp 
partial-wave analysis is somewhat smaller than 0.075. In the last line of Table | we quote 
f^o = 0.0745(6). 

In order to demonstrate that it is not only the data at low energies that pin down 
the coupling constant, Table |V1J| also contains the results for the ppir coupling constant 
obtained from a number of analyses in different energy ranges. It can be seen from this Table 
that the data at energies higher than 10 MeV or 30 MeV favor a low coupling constant as 
well. Similar results are found if we restrict the energy range at the high end, by doing an 
analysis of the data up to, say, 280 MeV. These findings once more underline our claim that 
the database as a whole contributes to a low value for the pion-nucleon coupling constant, 
and not some particular experiment (s) or the data in a restricted energy bin. 



V. pp ANALYZING-POWER DATA 
AROUND 10 MeV 

Let us next turn to our investigations of specific experiments that are discussed in con- 
nection with the pion-nucleon coupling constant. In Ref. |19|], for instance, it is stated that 
the prime reason for a low pprc coupling constant was our analysis of pp analyzing-power 



data around 10 MeV. In Ref. the same statement can be found in a different form, 
where it is said that the 3 P phase shifts around 10 MeV are very important in the deter- 
mination. However, these arguments only reflect our statements based on a preliminary pp 
analysis by our group 0. Apparently, these critics have overlooked our amendment to these 
statements as discussed in our paper on the completed 0-350 MeV pp analysis |l|], where we 
incorporated many theoretical improvements and included much more experimental data. 
In this latter paper it is explicitly stated that the 3 P waves are not especially important 
in the determination of f ppn0 and that it is not possible to pinpoint some specific type of 
observables as particularly constraining. Concerning the analyzing-power data around 10 



MeV, the 15 Wisconsin data points at 9.85 MeV p5] have in our latest multienergy pp anal- 
ysis Xmin = 16. If we fix the ppn coupling constant at 0.079 and refit, x m in on these data 
increases to 31. If we leave out this group (lowering the number of degrees of freedom with 
15 to 1598), x m in drops about 16 from 1786 to 1770 and the value for the coupling constant 
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becomes /p p7r o = 0.0751(6). This clearly shows that these data are not alone responsible for 
the low value of the coupling constant, although this group clearly favors a low coupling 
constant. We stress once more, however, that this latter conclusion is only justified when 
reached in a multienergy partial-wave analysis using the data as a whole. We did already 
demonstrate that the data above 30 MeV give for the coupling constant fp pn o = 0.0743(6), 
so the analyzing-power data around 10 MeV are absolutely not crucial to a low value for 
f 2 

A criticism of our pp analysis in this context is the fact that we do not include in our 
database another group of analyzing-power data around 10 MeV (and 25 MeV) taken by 
the Erlangen group of Kretschmer et al. |26|,|27| . We do not do this because it is our policy 
not to include data that have not been published in a regular physics journal. Moreover, 
in this specific case we have committed ourselves to not publishing any analysis of these 
specific data prior to their publication by the Erlangen group. Of course we are well aware 
of the existence of these data and we did analyze them. We can state that we find no reason 
whatsoever to modify any of our conclusions regarding the ppn° coupling constant. It is 
definitely not true that it is crucial which one of these two data sets (the Wisconsin or the 
Erlangen set) around 10 MeV is included, as is concluded in Ref. |19| from a study with 
meson-exchange potential models that have Xmin/^dat ~ 2. 



VI. np BACKWARD DIFFERENTIAL 
CROSS SECTIONS 

A point of discussion regarding the np analysis is the normalizations of the np differential 
cross sections and their relation to our determination of the np7r ± coupling constant (see, 
for instance, Ref. [0]). It is common folklore that the npn^ coupling constant is determined 
mainly by the peak present in backward np differential cross sections. It was suggested by 
Chew |28|] as early as 1958 that this is a good place to extract the pion-nucleon coupling 



constant. If this is true, then a very important group of data should be the Los Alamos set 



of backward cross sections measured by Bonner et al. |[29|| . However, we have seen already 
from Table [V] that all observables, and not the differential cross sections in particular, favor 
a low coupling constant. Our results and conclusions regarding the data from Bonner et 
al. are the following. The way we handle the normalization of a group of data and its 



uncertainty is explained in detail in Ref. p0[ . One group of 42 cross sections at 194.5 MeV 
is rejected completely, as well as 1 data point at 344.3 MeV. For the remaining 607 backward 
np cross sections at 10 different energies we find Xmm = 630 for / 2 = 0.075 and xLin = 655 
for / 2 = 0.079, so Ax^m = 25. Of these 10 groups, 7 groups have a floated normalization 
and 3 groups have a finite normalization error of 4%, where the sensitivity of these last 3 
groups (242 data points) to / 2 is rather small. So we see that it is mainly the shape of the 
cross section and not the normalization that makes that these Los Alamos data favor a low 



coupling constant. In Table [VIII| we give the results obtained for the charged-pion coupling 
constants for the 10 individual groups by fitting a parabola through 4 values of xLm f° r 4 
different coupling constants. We see that most groups favor a low coupling constant, but 
the errors are rather large. 

We have also tabulated the norm of these groups as determined in the multienergy 
analysis, once again following the x 2 -rise-by-one rule using the full error matrix. It can be 



S 



seen that our multienergy solution pins down these normalizations with very small errors, of 
the order of 0.5%, which is much smaller than the 4% error quoted by the experimentalists. 
There is essentially no difference between the groups with a floated normalization and the 
groups with a finite normalization error. 

Our conclusion is that the relevance of the backward cross sections for determining the 
charged-pion coupling constant is more limited than is generally assumed. There are other 
experiments that are much more constraining for the coupling constant. Here we mention the 



12 analyzing-power data at 10.03 MeV taken by Holslin et al. [pO] , the 16 spin correlations 
measured by Bandyopadhyay et al. [31] at 220 MeV, and the 19 spin correlations taken by 
the same group at 325 MeV. Again, these statements apply to an analysis of the groups 
within multienergy partial-wave analyses of the complete database, and do not follow from 
studies of the individual experiments. However, we want to stress once more that our 
low value of the charged-pion coupling constant is not only due to these accurate analyzing- 
power and spin-correlation data. The analysis without these data still yields / c 2 = 0.0750(4), 
demonstrating that also the other data favor a low, but slightly less accurate, value. 



VII. DETERMINATION FROM 
CHARGE-EXCHANGE DATA 

In this section, we add some remarks about our coupled-channels partial-wave analysis |§ 
of antiproton scattering data. In this case a neutral pion can be exchanged in elastic pp — > pp 
scattering and a charged pion in charge-exchange pp — > nn scattering. Until recently it was 
believed by probably everybody (including ourselves) that a partial-wave analysis of these 
reactions was out of the question. We find it gratifying that the methods used in the partial- 
wave analyses of pp and np scattering data could be extended to the case of the antiproton 
elastic and charge-exchange scattering. The value for the npn^ 1 coupling constant / c 2 = 
0.0751(17) found in our analyses of charge-exchange data is in nice agreement with the 
values found in the analyses of NN data. In fact, if it is possible to measure the differential 
cross section for pp — > nn with the accuracy stated by Bradamante (private communication, 



see also Ref. |j3~2|), the charge-exchange reaction will be an even more competitive place to 
study the isovector-meson coupling constants. 

We want to stress that the still popular (but now rather outdated) few-parameter optical- 
potential models can in no way be compared to a sophisticated partial-wave analysis. In 
the first approach at best a crude qualitative description of a limited number of data is 
possible. No Xmm * s ever presented. It is easy for us to construct a similar optical-potential 
model, by supplementing the C-parity-transformed Nijmegen potential |2TJ by an imaginary 
potential containing 2 free parameters. We then find at best Xmm ~ 10 9 f° r a database of 
3309 observables. This should be compared to xLin = 3592.5 reached in our multienergy 
partial-wave analysis on the same set of data. Of course, our two-parameter model is as 
bad or as good as any other few-parameter optical-potential model. For instance, for the 



1968 prototype Bryan-Phillips model [^] we find a Xmm which is even much larger. One 
can never hope to describe all pp scattering data with just 2 or 3 free parameters, when one 
needs already about 20 free parameters to fit the pp data. In a single-energy pp analysis in 
principle 8 times as many phase-shift parameters are required compared to a pp analysis. 
One should keep in mind that in the past these optical-potential models were never intended 
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for a quantitative comparison to the data. The 1980 Dover- Richard model ||34j| , for instance, 
served an excellent purpose in examining what could be expected qualitatively when LEAR 
would come into operation in 1983. At present, however, this approach seems hardly justified 
anymore. In our opinion, these naive models that do not fit the presently available data 
at all are completely inadequate to address in a reliable manner issues like the value of 
the pion-nucleon coupling constant, as was attempted very recently in Ref. |19 |. This can 



be seen, for instance, from the bad fit in Ref. [T9|] to very recent accurate charge-exchange 
analyzing-power data from LEAR ||35|| . After almost 10 years of data-taking at LEAR, it 
unfortunately still is a common practice to compare the data to the "predictions" of these 
museum models and then draw strong conclusions about the physics behind these models 
from such a comparison. 



VIII. CONCLUSIONS 

To summarize, we firmly believe that the value of the pion-nucleon constant found in 
the Nijmegen partial-wave analyses of pp, np, and pp scattering data is essentially correct 
and free of significant systematic errors. An excellent Xmm * s reached in all cases, reflecting 
both the statistical consistency of the data sets and the quality of the analyses. The specific 
method of analysis allows the extraction of the coupling constant at the pion pole from 
the asymptotic pion-exchange potential and ensures a clean separation from short-range 
form-factor effects and heavy- or multi-meson-exchange forces. We stress that in all cases 
we have also determined the mass of the exchanged pion and always found agreement with 
the experimental values. For instance, in the combined analysis of pp and np data it 
was found that m n o = 135.6(1.3) MeV and m n ± = 139.4(1.0) MeV. This success is a very 
strong argument against the presence of significant systematic errors, such as form-factor 
effects. Furthermore, there is consistency in the results from different analyses as well as 
agreement with the value f% = 0.0735(15) found by Arndt and coworkers in their latest 
VPI&SU analysis of n^p scattering data P,|36|,|3T|1 . 

We pointed out before || that the Goldberger-Treiman relation |38] also favors a low 



pion-nucleon coupling constant. Our present results indicate a need for reconsideration 
of calculations on the so-called Goldberger-Treiman discrepancy (see, e.g., Ref. |16|]). Re- 
cently, Workman, Arndt, and Pavan |39|] showed that the Goldberger-Miyazawa-Oehme sum 



rule [[K]] provides rather model-independent evidence for a low coupling constant as well. 
In a recent study of the deuteron properties |T7| it was shown that a low value for 



the pion-nucleon coupling constant implies that the value for n p = /awp/ 9nnp = 6.6 as 
determined in the Karlsruhe-Helsinki analyses of n^p scattering data |]TU[ must be wrong. It 
was further shown that the preferred value for k p is in agreement with the value k p = 4.2 as 
found in the 1978 Nijmegen soft-core nucleon-nucleon potential [21] and with the value k p 



= 3.7 which follows from vector-meson dominance of nucleon electromagnetic form factors. 

We feel that especially the value for the ppn° coupling constant as determined in the pp 
analysis is compelling evidence for a low pion-nucleon coupling constant. With the excep- 
tion of the 1 Sq wave, the ppn coupling constant can be determined from all parametrized 
partial waves, all values being consistent. Given the ppn° coupling constant, it seems to us 
that claims for a high npn^ 1 coupling constant are untenable, since in that case not only 
three independent recent determinations of this coupling constant must be wrong, but one 
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also has to cope with a large breaking of charge-independence, which is theoretically very 
difficult to accommodate ||. Therefore, we strongly recommend that in future work on 
nucleon-nucleon scattering the value f% NlT = 0.0745 at the pion pole is taken as a starting 
point from which further consequences can be discussed. 
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TABLES 



TABLE I. Recent determinations of the pion-nucleon coupling constants from disper- 
sion-relation (DR) analyses of pion-nucleon scattering data and from partial-wave analyses (PWA) 
of nucleon-nucleon and antinucleon-nucleon scattering data. 



Group 


Year 


Method 


10 3 / 2 o 


10 3 / 2 


Karlsruhe-Helsinki fl2]| 


pre-1983 


ir^~p DR 




79(1) 


Nijmegen [jl|] 


1987-1990 


pp PWA 


74.9(0.7) 




VPI&SU § 


1990 


it^p DR 




73.5(1.5) 


Nijmegen 0| 


1991 


combined NN PWA 


75.1(0.6) 


74.1(0.5) 


Nijmegen || 


1991 


pp PWA 




75.1(1.7) 


this work 


1992 


pp and np PWA 


74.5(0.6) 


74.8(0.3) 



TABLE II. Comparison between the moments of the x 2 -P r obability distribution expected from 
theory and those determined in our partial-wave analysis (PWA) of pp data. Tabulated are 
(x 2 nin )/iVdf and the central moments /i n for n = 2,3,4. 





theory 




PWA 


(X 2 min )/N di 


1.000±0.035 




1.108 




1.81±0.12 




1.83 


M3 


5.55±0.74 




5.40 


M4 


29.8±4.5 




27.6 


TABLE III. The ppir coupling 


constant as a function of the cutoff mass in 


the exponential 


form factor. The number of degrees 


of freedom is 1612. 






A (MeV) 500.0 


750.0 1000.0 


1250.0 


oo 


103 O 75.2(0.5) 


75.0(0.5) 75.0(0.5) 


75.0(0.5) 


75.0(0.5) 


Xlin 1786.3 


1786.4 1786.4 


1786.4 


1786.4 
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TABLE IV. x 2 results for 4 values of the ppir coupling constant for the different types of 
observables in the analysis of pp scattering data. The numbers in the two last columns are obtained 



Irv fitting 


a parabola to the numbers 


in the four preceding 


roll limns 


The number of degrees of 


freedom is 


1613 for each value of f 2 

J ppn u 












type 


iVdat 


10 3 O = 73 


75 


77 


79 


X (mm J 


10 3 f 2 o(min) 


da /dfi 


821 


838.7 


825.7 


823.0 


830.8 


822.7 


76.5(0.9) 


Ay 


558 


585.2 


580.0 


585.7 


602.2 


580.0 


75.0(0.9) 


a ■■ r 


66 


52.8 


55.5 


60.0 


66.3 


51.9 


71.0(2.1) 


D,D t 


97 


105.2 


107.5 


112.0 


118.9 


104.9 


72.0(1.9) 


R,R',A,A' 


209 


194.5 


193.1 


194.9 


199.8 


193.1 


74.9(1.6) 


rest 


36 


24.8 


24.7 


24.6 


24.5 






all 


1787 


1801.2 


1786.4 


1800.2 


1842.4 


1786.4 


75.0(0.5) 



TABLE V. x 2 results for 4 values of the npit^ coupling constant for the different types of 
observables in the analysis of np scattering data. The numbers in the two last columns are obtained 
by fitting a parabola to the numbers in the four preceding columns. The NNt: coupling constants 
are taken to be 0.075. The number of degrees of freedom is 2331 for each value of f 2 . 



type 


iVdat 


10 3 f 2 = 73 


75 


77 


79 


X 2 (min) 


10 3 / c 2 (min 


(7 t ot,Ao"L,A(7T 


252 


232.9 


229.7 


232.4 


242.4 


229.5 


75.1(1.1) 


da/dn 


1350 


1379.0 


1364.2 


1367.9 


1391.8 


1363.2 


75.6(0.6) 


Ay 


738 


737.7 


720.3 


745.9 


830.4 


717.8 


74.8(0.4) 


A A 


86 


77.2 


72.6 


91.2 


136.0 


71.2 


74.4(0.6) 


D t 


43 


42.8 


39.8 


42.0 


51.6 


39.5 


75.1(1.1) 




43 


54.6 


58.6 


68.5 


88.4 


54.7 


73.1(1.0) 


all 


2512 


2524.3 


2485.3 


2547.9 


2740.7 


2480.4 


74.8(0.3) 



TABLE VI. Results for the pion-nucleon coupling constant introduced separately in each 
parametrized partial wave in the analysis of pp scattering data. For fitting the coupling con- 
stant in non-S waves the coupling in the 1 S'o wave is fixed at 0.075. The number of degrees of 
freedom is 1611 in each case. 



partial wave 


10 3 f 2 (wave) 


10 3 / 2 (rest) 


X 2 - 


'So 


79.7(1.9) 


74.5(0.6) 


1779.8 


l D 2 


74.6(0.8) 


74.6(0.6) 


1786.0 




74.6(2.1) 


74.6(0.6) 


1786.0 




72.7(1.7) 


74.8(0.6) 


1784.6 


"Pi 


74.9(0.7) 


74.3(0.8) 


1785.6 


3 P 2 - 3 F 2 


74.7(0.8) 


74.6(0.6) 


1786.0 


3 F 3 


73.3(1.3) 


74.8(0.6) 


1784.8 


3 F 4 - 3 F 4 


75.1(0.9) 


74.5(0.6) 


1785.6 
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TABLE VII. Values for the ppir coupling constant determined in a partial-wave analysis pp 



scattering data within different energy ranges in MeV. 

flab range 10 3 f 2 ppn o Xmin ^min 

/N di 

0-350 1612 75.0(0.5) 1786.4 1.108 

3-350 1435 74.3(0.6) 1596.6 1.113 

10-350 1312 73.7(0.7) 1488.8 1.135 

30-350 1237 74.2(0.8) 1397.6 1.130 

0-280 1243 75.5(0.6) 1389.3 1.118 

3-280 1066 74.5(0.7) 1189.7 1.116 



TABLE VIII. The charged-pion coupling constant f£ determined from the backward np cross 



sections of Bonner et al. 2£] in the np partial- wave analysis. The numbers are obtained by fitting 
a parabola through the Xmin results for 4 different coupling constants. For Ti a b = 265.8 MeV these 
4 numbers were consistent with a straight line. 



flat (MeV) 



at 



X 2 (min) 



norm 



10 3 / c 2 (min) 



162.0 


43 


60.0 


1.092(7) 


69.9(3.0) 


177.9 


44 


44.0 


1.083(7) 


70.2(3.1) 


211.5 


43 


31.0 


1.063(7) 


72.8(3.3) 


229.1 


49 


62.3 


1.058(7) 


69.5(3.6) 


247.2 


53 


38.5 


1.042(7) 


69.7(9.3) 


265.8 


63 




1.028(6) 




284.8 


73 


79.7 


1.052(5) 


75.3(3.5) 


304.2 


80 


79.9 


1.003(4) 


74.6(3.4) 


324.1 


82 


91.7 


1.057(5) 


78.0(4.3) 


344.3 


80 


74.6 


1.035(5) 


74.8(3.9) 
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